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I . INTRODUCTION 

In the design of hypersonic aircraft, the relatively high heat- 
ing loads encountered make it difficult to provide an engine 
with a nozzle of sufficient length to achieve complete internal 
expansion. In this case, optimum performance of the system 
dictates cutting several walls of the nozzle short in order to 
achieve a savings in weight and making use of the vehicle under- 
surface as the final nozzle expansion surface. Two such hyper- 
sonic vehicles are represented schematically in Figure (1). In 
each case, the vehicle is propelled by a number of engines 
which all discharge through relatively short nozzles into a 
final single common nozzle formed to a large extent by the ve- 
hicle undersurface. 

Since the upper wall of the common nozzle is formed by the lower 
surface of the vehicle its length is fixed by vehicle design 
considerations rather than by considerations of nozzle optimiza- 
tion. The optimum length for the nozzle cannot be defined with- 
out a series of calculations since its value is determined by a 
compromise between increased efficiency and increased structural 
weight. In order to perform the nozzle optimization, it is nec- 
essary to carry out a series of aerodynamic analyses which fur- 
nish a family of designs and their corresponding performance 
and to determine from these the relation between length, and 
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therefore weight, and perform ance . It is the purpose of this 
volume to outline such a procedure. 

In practical applications, the flow leaving the burners of 
ramjet and scranjet engines is far from uniform with large 
variations of stagnation conditions and flow prooerties exist- 
ing across the burner exit plane. Moreover, the flow at the 
nozzle exit is also nonuniform, since the exit area is smaller 
than that which is required for full expansion and the nozzle 
length is necessarily shorter than that required for comDlete 
wave cancellation. As a result of this overall nonunifornity, 
it is not possible to develop exact criteria for the design 
of an optimum nozzle. Fortunately, the sensitivity of nozzle 
performance to minor perturbations of exit Mach number and 
flow inclination is quite small. Therefore, satisfactory 
nozzle designs may be obtained by means of an iteration pro- 
cedure where initial values of several of the design parameters 
are assumed from simplified approximate considerations and are 
then varied systematical ly in order to determine the value of 
each parameter which optimizes nozzle performance. 

This procedure can be carried out using the three dimensional 
nozzle analysis described in Volume I of this report. The 
geometric design parameters are initially selected from euasi- 
two dimensional considerations which will be described in the 
following sections and the performance of the selected con- 
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figuration is determined using the computer programs described 
in Volume II of this report. Then, improvements in nozzle 
performance are obtained by perturbing the initially assumed 
values of the geometric parameters on the basis of the results 
of the computer analysis, this procedure being followed until 
no further performance improvements are obtained. The results 
of this study can then be used to evaluate the trade-off be- 
tween increased nozzle efficiency and increased nozzle weight. 
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II. DESCRIPTION OF NOZZLE REQUIREMENTS 


The geometry of a general multiple nozzle assembly is outlined 
schematically in Figure (2). 



FIGURE 2. BOTTOM AND REAR VIEW OF GENERAL MULTIPLE NOZZLE 
CONFIGURATION 

The walls of the nozzles expand simultaneously in two directions 
approximately normal to one another. Since the burner design 
fixes the thickness (t) of the wails which divide the individual 
nozzles, it therefore fixes the degree of lateral expansion 
which will be produced when these walls are removed. 
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Because the principal expansions occur on surfaces which are 
approximately normal to one another, the net change in Mach 
number produced can be accurately predicted by two dimensional 
theory using the sum of the turning produced by each individ- 
ual wave as the total turning angle. The shapes of the waves, 
of course, cannot be predicted in this way and must be deter- 
mined by the nozzle analysis program. 

The surfaces BC and EF are planes whose orientation is known. 

The nozzle analysis does not require these p.lanes to be paral- 
lel but their inclination must be fixed. As a first step, it 
is assumed that the wall HI is also given and planar until the 
station at which the solid surfaces AB, DE and 3H are termi- 
nated. Other geometric boundary conditions could be imposed 
on the analysis, but these would only produce losses in thrust 
since lateral external expansion would take place before the 
flow was sufficiently expanded in the nozzle and would require 
that the end nodules be calculated separately. Therefore, 
lateral expansions are considered only downstream of point H 
(surface HI is present, but I can coincide with H). 

Since the nozzles are three dimensional, they expand in vertical 
as well as lateral planes. The general geometry of the nozzle 
in the vertical plane is outlined in Figure (3). 
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A* 

A 

FIGURE 3. SIDE VIEW OF GENERAL MULTIPLE NOZZLE CONFIGURATION 

The multiple primary nozzles terminate along a curved line BB‘. 
The cowl surface, AC, Is far shorter than the undersurface of 
the fuselage, A'J. The side walls of the nozzle GI, terminate 
along the line CC which is inclined so as to form a super- 
sonic trailing edge. The shape of the vehicle undersurface, 

C'J, is assumed to be given. 

In practical aircraft designs, the length of the cowl surface 
AC is far shorter than that required to provide the possibility 
of cancelling the expansion waves eminating from A', while the 
vehicle length A'J is sufficiently large so that at least a por- 
tion of the waves produced at A can be cancelled. In general. 
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the surfaces AB and A'B 1 are not sufficiently long to feel 
the effects of the expansion fans from the corners A and A'. 
Therefore, wall curvature upstream of the line BB' would only 
produce premature compression waves and reduce nozzle effi- 
ciency. As a result, these short surfaces are taken as 
straight walls and are completely characterized by their in- 
clination to the vehicle centerline. 

Summarizing, subject to the constraints of a given length and 
vertical rise for the vehicle undersurface A'O, and a given 
thickness, t, of the sidewalls of the multiple nozzles, the 

I 

geometry of the system can be defined once the inclinations 

of the surfaces A *B * and AB and the sidewall shapes (AB, 

DE, etc.) have been specified. The cowl length, AC, can then 
be reduced and the resulting nozzle losses computed for a 
comparison with the weight savings effected. 

It should be noted here, that since the waves from A and A' are 
not completely cancelled at the opposite walls, the flow pro- 
perties along the line CJ are far from uniform. In this region, 
those waves produced along AA' are still present but are suf- 
ficiently weak to allow two dimensional criteria to be applied. 
This approximation will be used to select the inclinations of 
the surfaces AB and A'B'. 
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III. SELECTION CRITERIA FOR -HE INITIAL NOZZLE GEOMETRY 

Since the actual area available for the nozzle exit is usually 
insufficient to allow either complete expansion or complete 
wave cancellation, the nozzle exit flow will be nonuniform and 
at an average static pressure somewhat higher than ambient. 

The general wave pattern which will exist in a vertical refer- 
ence plane through the nozzle is indicated in Figure (4) below 



FIGURE 4. 


WAVE PATTERN IN VERTICAL NOZZLE REFERENCE PLANE 
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The difference in static pressure between the wall at 0 and 
the central streamtube is determined by the initial flow Mach 
number, the geometry and the degree of turning at A, the com- 
bination of which determines how much of the expansion fan 
generated at A reaches the wall A'J, and by the degree of wave 
cancellation along the wall. Similarly, the pressure differ- 
ence between the cowl exit point C and the central streamtube 
is determined by the degree of turning at A* and the amount 
of wave cancellation along the cowl AC. In practical designs 
the length of the cowl is usually so short as to make the 
degree of cancellation along its surface small enough to be 
neglected. 

Two important nozzle properties are determined by overall ve- 
hicle design considerations: the total nozzle exit area is 
fixed by vehicle geometric constraints and the angle of the 
thrust vector with respect to the free stream is determined 
by lift and pitching moment requirements. In addition, the 
initial flow inclination 6 is a function of inlet and combustor 
design considerations and is assumed to be a nrespecified para- 
meter. 

The first step in the nozzle design orocess is the definition 
o-f the entrance plane flow conditions. Since the nature of 
ramjet and scramjet engines is such that the combustor exit 
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flow is usually highly nonuniform in both stagnation pressure 
and enthalpy, a simple treatment of the nozzle flow field re- 
quires that the entrance conditions be idealized somewhat. 
This is accomplished in the present case by dividing the noz- 
zle entrance plane into regions as indicated in Figure (5). 

In each of these regions, an engineering approximation to the 
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FIGURE 5. DIVISION OF THE NOZZLE ENTRANCE PLANE 

total enthalpy may be obtained by dividing the area into a 
number of subregions of nearly constant properties and mass 
averaging the total enthalpies of each subregion. To obtain 
the^average total pressure, the procedure is also straight- 
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forward but slightly more involved numerically. In each sub- 
region, i, we must first determine the flach number M-j, the mass 
flow m-,* and the total temperature T T _.. Then, the mass flow 
function N ^ is calculated for each subregion using the defini- 
tion 

-1 

Nj = M ^ ( 1 + M^) ( 1 + y M ? ) (1) 


The average mass flow function for the entire program is then 
given by Reference (1) to be 


N 


n 

Z 

i = 1 


/m -i Z 
i = l 




n 

i*l T. 


/mi Tj. 



( 2 ) 


Knowing fT, the value of TT is found by a solution of Equation (1). 
Assuming a constant value of static pressure, TT, across the noz- 
zle entrance plane the average total pressure for the region is 
given by 


p r = p (1 


Y- 1 

+ 2 M ) 


_X_ 

1-Y 


( 3 ) 


Having defined the necessary average properties, each of the 
regions, is now expanded by means of a one dimensional analysis 
to a pressure no less than - the local ambient value. The values 
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of streamtu.be area, A, and gas velocity, V, are determined 
as functionsof static pressure through the expansion as is 
the increase in two dimensional Prandtl -Meyer turning angle, 
v, corresponding to the instantaneous value of Mach number. 

Since the purpose of this analysis is simply to determine 
starting values for the nozzle geometric parameters used in 
the iteration procedure described in the previous section, a 
high degree of accuracy is not reouired. As a result, several 
assumptions have been and will be introduced in order to 
simplify the calculation. At this point, another simplifica- 
tion may be properly introduced: the value of y may be as- 
sumed constant within prescribed ranges of temperature, jump- 
ing disconti nuously at the boundary from its value in one 
range to its value in the next. By so doing, existing tables 
of one dimensional flow properties may be utilized for this 
portion of the calculation without a significant sacrifice in 
accuracy. The result of this computation is shown schematical- 
ly in Fi gure (6 ) . 



The total increase in Prandtl -Meyer turning angle experienced 
by each of the three streams is simply 

Av a = Av^ + Av A i + t 

Av b = Av l + Av a + Av a , 

Av c = Av^ + Av^ 

where Av^ and Av^, are the degrees of turning produced by the 
corners A and A', Av L is the degree of lateral expansion* and 
t is the degree of turning produced by the waves from A' which 
reach the vehicle undersurface as indicated in Figure (4). 

Denoting the i ncl inati on of the thrust vectors for the three 
regions with respect to the free stream direction 1 as o a , a b 
and a c> geometry dictates that 

a a = 6 + Av^ ( - t 

= S + AV ^ , - Av A 

® c * 6 • fv /> 

*Av»l '"i's calculated by takYng' the change in Prandtl -Meyer angle 
corresponding to the expansion from the initial entrance ilach 
number produced by the lateral area ratio (d+t)/d. 
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Since the inclination, or of the resultant thrust vector is 
a given quantity, dictated by aircraft design considerations, 
we require that 


( a a" a R ) ( m a v a + Pa A a ) + ( a b _a R) ( m b v b + Pb A b) + 
(Oc-'r) ("cVl’cV " 0 


( 10 ) 


Finally, since the total area. A, at the nozzle exit plane is 
fixed, we have the condition 


A, + Au + A„ = A 


( 11 ) 


Mow, several values of the exit plane static nressure in strean- 

tube b are selected, either equal to or greater than the local 

ambient static pressure. For each pressure selected, the 

total flow deviation Av b ,, the velocity V b , and the streamtube 

area A b , are given by the one dimensional calculation represented 

in Figure (6). Using these values, Equations (4) through (11) 

are used to provide a numerical solution for the corresponding 

values of Av f , Av^ ( and t. The operation is repeated for 

several values of the pressure n' and the optimum set of values 

b 

is selected for the initial guess as to nozzle geometry. 


As stated earlier, since the length of the multiple nozzle 
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section is relatively short very few, if any, of the waves 
generated at the corners A and A' will reach the cowl or 
vehicle surfaces upstream of the sidewall trailing edges CC' . 
As a result, these upper and lower surfaces should be kept 
linear within the confines of the multiple nozzles since pre- 
mature curvature of these walls would produce envelope shocks 
and increase losses. However, the nozzle sidewalls must be 
reduced in thickness from their initial value of t to zero 
within the distance from A to C. Since, directional con- 
tinuity at the trailing edge requires the formation of shock 



FIGURE 7. EFFECT OF STRAIGHT SIDEWALLS 

waves of strength 6, restriction of the trailing shock losses 
to acceptable limits requires that the trailing edge angle 
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be quite small. Thus, if the sidewalls were to be designed 
as straight surfaces and the shock losses were to be kept 
small, the sidewall length would become excessive. To avoid 
this problem, curvature of the sidewalls is indicated and 
the degree of this curvature becomes another design parameter. 

Consider the lateral wave pattern depicted schematically in 
Figure (8). By turning the sidewalls AB and PE gradually 



FIGURE 8. WAVE PATTERN* WITH CURVED SIDEWALLS 
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back toward the axial direction, a series of compression 
waves are generated which coalesce downstream to form an 
envelope shock. Although the eventual strength of this wave, 
when fully formed, cannot exceed the strength <5 of the wave form 
ed by a straight wall design, we wish to design the wall to 
delay coalescence for the maximum possible distance. In 
this regard, the vertical expansion waves eminating from 
the corners A and A' are quite helpful, since they diffract 
the waves and tend to move the region of coalescence down- 
stream. 

Therefore, the design process should proceed along the follow- 
ing lines. For each of the three vertical regions, a, b and 
c defined earlier, the Mach numbers both upstream and down- 
stream of the expansion fan from corners A and A 1 are known. 

In each of these regions, the sidewalls are designed using 
two dimensional characteristics techniques as described in 
Reference (2). For the two regions a and c bounded by the 
cowl and vehicle surfaces, the sidewalls are designed using 
the Mach number and flow deflection downstream of the corner 

expansion as the initial conditions: in region b, the central 

core, the combustor exit Mach number and direction are more 
appropriate. Once the three profiles are defined, the over- 
all geometry may be defined by a linear variation with verti- 
cal rise from profile c to profile b to profile a. It should 



TR 166 - Vol . Ill 


Page 19 


be noted that it is not necessary to utilize the averaged 
properties in each region for the characteri si tcs computa- 
tion since the methods presented in Reference (2) allow 
handling discontinuities in entropy and enthalpy across the 
f 1 o w field. 


As a guide to the proper range of values for the initial 
sidewall expansion angle e, consider the total amount of 
lateral expansion which must be produced as defined by the 
area ratio (d+2t)/d. Starting from the initial Mach number, 
this area ratio corresponds to a specific increase in the 
Prandtl -Meyer function Av L . This total change in v is 
exactly twice the initial expansion angle reouired if uniform 
shock free flow were to be obtained in the lateral nlane. 

This value will generally be fairly small and will lead to 
a rather long sidewall. Initial choices of larger lateral 
expansion angles, c, will nroduce shorter walls and increas- 
ing wave coalescence. Summarizing 


e > 



( 12 ) 


As noted earlier, the three dimensional effects in the nozzle 
diffract the sidewall compression waves and delay formation 
of the envelone shock. Therefore, the sidewall geometry de- 
signed using two dimensional technioues will be longer than 
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actually required. As a consequence, once the flow field for 
the selected initial geometry is analyzed using the program 
described in Volume II of this report, the sidewall geometry 
may be perturbed in the direction of larger angles e and cor- 
respondingly shorter overall lengths 



